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INTRODUCTION 


1  . 

1.1  OVERVIEW 

The  ability  to  conduct  long  term  experiments  in  space  has  necessitated 
understanding  the  local  environment  of  the  measuring  platform  in  order  to 
be  able  to  sort  out  the  perturbations  generated  by  the  measuring  platform 
on  the  phenomena  being  measured.  A  number  of  attempts  have  been  made  to 
model,  in  a  piecemeal  fashion,  the  effect  of  contamination  on  sensors  flown 
on  the  Space  Shuttle.  Many  of  the  earlier  attempts  have  treated 
contamination  as  being  due  solely,  or  primarily,  to  the  exhaust  products  of 
thruster  motors. Of  course,  spacecraft  contamination  is  indeed 
related  to  the  flow  fields  of  the  exhaust  products  of  spacecraft  engines, 
as  has  been  demonstrated  in  a  Monte  Carlo  code  developed  for  that  specific 
purpose  of  understanding  the  backflow  of  exhaust  products  into  the  shuttle 
bay. This  report  describes  calculations  performed  with  SOCRATES 
(Spacecraf t/Orbiter  Contamination  Representation  Accounting  for  Transiently 
Emitted  Species)  which  has  been  developed  with  the  idea  of  addressing  in  a 
comprehensive  manner  the  problem  of  contamination  on  spacecraft.  SOCRATES 
builds  on  a  plume  code  developed  for  a  different  application^)  and  tries 
to  assess  the  flow  fields  of  contaminants  in  the  low  earth  orbit 
environment  taking  into  account  not  only  scattering  but  reactive 
collisions.  The  importance  of  this  effort  arises  from  the  fact  that  at 
orbital  velocity  in  low  earth  orbit  (7,7  km  s'*)  spacecraft  surfaces  as 
well  as  the  gaseous  envelope  of  the  spacecraft  are  constantly  bombarded  by 
atomic  oxygen.  The  velocity  distribution  of  the  0  atoms  is  typical  of 
ionospheric  temperature,  measured  as  750  ±  50  K  on  one  occasion.  (6)  At 
orbital  velocity  the  energy  of  0  is  4.9  eV,  so  the  potential  for  undergoing 
endothermic  reactions  is  great.  The  timeliness  of  code  development  is 
particularly  propitious  because  of  the  availability  of  data  from  some  Space 
Shuttle  flights  where  diagnostic  instruments,  e.g.  mass  spectrometers (7-)^) 
and  plasma  diagnostics^ )^)  were  flown.  It  has  been  clearly  shown  in  these 
experiments  that  gases  from  some  of  the  exhaust  jets  are  reflected  into  the 


exhaust  jets  are  reflected  into  the  bay  of  the  Space  Shuttle,  the  amount 
being  dependent  on  the  orientation  of  the  VRCS  with  respect  to  the  velocity 
vector  and  with  respect  to  the  shuttle  surf aces. 

The  eventual  sin  in  the  development  of  SOCRATES  is  the  inclusion  of 
gas  phase  neutral-neutral  reactions,  neutral-surface  reactions,  and 
neutral-ion  reactions  and  to  validate  the  code  against  the  measurements 
reported  above  as  well  as  against  those  which  measured  surface 
interactions.  This  report  presents  calculations  performed  with  the 
first  module  of  SOCRATES,  which  is  the  Monte  Carlo  module  for  far  field 
neutral  gas  interactions. 

1 . 2  METHODOLOGY 

The  direct  simulation  Monte  Carlo  method  involves  storing  a  discrete 
number  of  molecules  (via  their  velocities,  positions,  and  other  pertinent 
information)  in  a  computer.  The  solution  region  is  broken  up  into  a  number 
of  separate  cells,  and  the  solution  is  stepped  forward  in  time  in  a  two 
stage  process.  First,  the  molecules  are  advanced  along  their  trajectories 
by  an  amount  appropriate  to  their  velocity  and  a  time  increment,  &tn.  In 
this  first  stage  some  molecules  will  leave  the  solution  region,  and  some 
will  be  introduced  as  determined  by  the  boundary  conditions  for  a 
particular  problem.  The  second  stage  is  to  simulate  collisions  in  each 
cell  appropriate  to  atm  so  that  collision  frequencies  are  properly 
simulated.  A  basic  hypothesis  of  the  method  is  that  if  the  time  step  is 
made  small  enough,  the  processes  of  translations  and  collisions  can  be 
uncoupled  in  this  manner. 

Periodically,  the  solution  is  sampled  by  accumulating  statistical  sums 
of  number  densities,  velocities  and  other  basic  properties.  The  solution 
is  run  repeatedly  until  statistical  deviations  are  reduced  to  a  desired 
limit,  and  then  physically  meaningful  output  quantities  are  computed  from 
the  statistical  sums.  The  number  of  molecules  represented  is  typically 
many  thousand  at  a  time,  which  is  vastly  fewer  than  the  number  occurring  in 
virtually  all  real  flows.  Hence,  the  construction  of  a  dynamically  similar 
flow  to  be  simulated  in  the  computer  is  an  essential  feature  of  the  method. 


In  order  to  successfully  apply  this  technique  to  the  contamination 
flow  fields  of  interest,  consideration  must  be  given  to  the  basic  character 
of  these  flows.  A  molecule  leaving  the  shuttle  can  be  expected  to 
experience  its  first  collision  with  an  atmospheric  species  at  a  distance, 
»c,  from  the  shuttle  which  depends  on  the  local  conditions.  Typically,  for 
unperturbed  space,  *c  may  be  on  the  order  of  a  few  hundred  meters  at  an 
altitude  of  200  km,  some  tens  of  kilometers  at  400  km,  and  hundreds  of 
kilometers  at  an  altitude  of  600  km.  The  solution  region  must  extend  to 
several  tc  if  the  scattering  process  is  to  be  accurately  represented,  which 
is  necessary  for  calculating  atmospheric  backscatter  or  radiative  emission 
resulting  from  these  collisions. 

A  fundamental  problem  arises  in  describing  accurately  the  interaction 
of  the  scattered  flow  with  the  shuttle  as  *c  becomes  large,  since  cell 
sizes  must  be  small  compared  to  shuttle  dimensions.  The  required  number  of 
cells  for  such  a  complete  solution  makes  the  complete  calculation 
intractable.  However,  advantage  can  be  taken  of  the  separation  of  length 
scales  to  achieve  the  same  result  with  much  less  computational  effort. 
This  is  done  by  separating  the  solution  into  an  "inner"  and  an  "outer" 
solution,  the  latter  being  performed  first. 

The  outer  solution  is  performed  on  a  length  scale  of  the  interaction 
of  the  contaminants  with  the  atmosphere,  i.e.,  several  *c.  On  this  length 
scale  the  dimensions  of  the  shuttle  are  unimportant,  and  the  sources  of 
contaminants  are  treated  as  emanating  from  the  origin.  There  are  no 
shuttle  surfaces  considered  in  the  outer  solution,  and  the  smallest  cells 
are  large  with  respect  to  shuttle  dimensions.  The  outer  solution  is 
generated  by  the  first  module  of  SOCRATES,  which  is  described  in  this 
paper . 

The  inner  solution  is  performed  by  the  second  module  of  SOCRATES, 
which  is  still  under  development.  The  inner  solution  comprises  the 
immediate  shuttle  vicinity,  and  takes  as  a  boundary  condition  the  results 
from  the  outer  solution  for  the  inwardly  directed  fluxes  of  the  species  in 
the  simulation.  The  inner  solution  handles  the  detailed  interaction  of  the 
contaminant  cloud  with  the  shuttle  vehicle,  as  well  as  direct  contamination 
from  sources  such  as  the  backflow  from  a  thruster.  In  the  inner  solution, 
all  sources  are  properly  located  on  the  vehicle. 
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Once  che  inner  and  outer  solution  have  been  finished,  then  ion-neutral 
interactions  will  be  included  in  the  code.  This  modification  will  require 
a  number  of  detailed  kinetic  data  as  well  as  inclusion  of  electric  fields 
into  the  model. 


2.  CASE  DEFINITIONS 

2.1  PHYSICAL  PARAMETERS 

Calculations  have  been  performed  for  the  interaction  between  the 
ambient  atmosphere  end  the  steady  state  exhaust  from  a  Space  Shuttle  VRCS 
(Vernier  Reaction  Control  System)  engine.  It  should  be  stressed  that  this 
choice  is  arbitrary  and  is  not  a  limit  of  SOCRATES.  A  suite  of  nine 
calculations  were  performed,  for  altitudes  of  200,  400,  and  600  km,  and 

angles  of  attack  between  the  oncoming  wind  ar.J  the  centerline  exhaust 
direction  of  0,  90,  and  180  degrees.  (A  180°  angle  of  attack  corresponds 
to  firing  the  VRCS  in  the  upwind  direction,  while  0°  corresponds  to  the 
downwind  direction. ) 

The  predicted  composition  of  the  exhaust  is  shown  in  Table  1  and  is 
based  on  thermodynamic  equilibrium  at  3000  K  (quoted  by  Pickett  et  al.^9^). 
The  constituents  of  the  exhaust  interact  with  the  ambient  atmosphere,  whose 
composition  is  given  in  Table  2. 

The  VRCS  engines  release  gases  at  a  velocity  which  is  typical  of 
chemical  rockets,  about  3  km/s.  Hence,  the  relative  velocity  between 
exhaust  and  atmospheric  species  varies  from  about  4.7  to  10.7  km/s, 
depending  on  the  angle  between  the  exhaust  and  atmospheric  velocity 
vectors.  For  atomic  oxygen  colliding  with  water,  for  example,  this 
corresponds  to  an  energy  range  of  about  1  to  5  eV.  Therefore,  a 
substantial  range  of  collision  energies  is  available  depending  on  the 

relative  velocity  between  the  exhaust  and  atmospheric  species. 

2.2  COORDINATE  SYSTEM 

A  Cartesian  coordinate  system  was  used  with  its  origin  at  the  shuttle 
and  the  +X  direction  fixed  as  the  direction  of  shuttle  motion.  This  means 
that  for  an  observer  fixed  with  respect  to  the  shuttle,  the  atmospheric 

wind  is  heading  in  the  -X  direction.  For  the  0°  angle  of  attack  cases  the 

centerline  exhaust  direction  is  also  in  the  -X  direction,  while  it  was 
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taken  to  be  in  the  +Z  direction  for  90°  and  +X  direction  for  180°  angles  of 
attack.  Hence,  the  X-Z  plane  is  a  symmetry  plane  for  all  cases,  while  the 
0°  and  180°  cases  are  also  axially  symmetric  with  respect  to  the  X  axis. 

2.3  FAR  FIELD  LINE-OF-SIGHT  INTEGRATIONS 

As  was  discussed  in  Section  1,  the  scale  lengths  for  the  interactions 
were  substantially  different  in  many  of  the  cases,  making  a  common  spatial 
scale  impractical.  Hence,  when  comparing  results  at  different  altitudes, 
the  reader  should  remember  that  the  presentation  length  scales  may  be  quite 
different.  In  many  cases  the  results  are  presented  as  integrated  along  the 
Y  direction,  which  shows  how  they  would  appear  to  a  distant  observer.  This 
is  mainly  done  in  order  to  reduce  the  dimensionality  of  the  presentation, 
while  retaining  substantial  physical  meaning.  The  presentation  of 
quantities  which  vary  as  a  function  of  three  spatial  dimensions  is  not 
straightforward  on  two  dimensional  paper. 

The  far  field  integrations  are  presented  in  two  forms:  gray  scale  and 
contour  plots.  The  gray  scale  plots  have  the  advantage  that  they  give  an 
immediate  visual  impression  of  the  relative  magnitude  of  the  spatial 
variation  of  a  quantity.  The  contour  plots  have  the  advantage  that  they 
can  be  more  quantitatively  read. 

2.4  NEAR  FIELD  LINE-OF-SIGHT  INTEGRATIONS 

Finally,  an  attempt  was  made  to  calculate  the  quantities  as  they  would 
be  seen  by  an  observer  located  within  the  shuttle  bay  as  a  function  of  look 
angle.  Such  calculations  can  only  be  done  accurately  after  the  inner 
solution  portion  of  SOCRATES  is  complete,  at  which  time  substantially  more 
detail  would  be  expected. 

When  a  particular  viewing  direction  is  considered  from  the  shuttle,  it 
is  represented  by  the  angles  a  and  0.  These  angles  are  illustrated  in 
Figure  1,  which  shows  an  arbitrary  viewing  direction  represented  by  the 
viewing  direction  unit  vector  fy.  a  is  the  angle  between  the  -OC  direction 
and  Iv,  and  0  is  the  angle  between  the  +Y  direction  and  the  projection  of 
iv  into  the  Y-Z  plane.  Hence,  iy  is  represented  in  terms  of  a  and  0  as 
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cos(a)Ix  +  sin(a)  cos  O)  iy  +  ain(a)cos(0)Iz 


(1) 


This  definition  was  chosen  since  for  the  axisymmetr ic  cases  there  is  no  0 
dependence  for  quantities  integrated  along  a  line  of  sight.  The  near  field 
results  will  also  be  shown  for  both  gray  scale  and  contour  plots. 


TABLE  1 

Shuttle  Engine  Combustion  Products1 


Species 

Mole 

Fraction 

Molecular 
Weight  (AMU) 

h2o 

0.33 

18 

n2 

0.31 

28 

h2 

0.  16 

2 

CO 

0. 13 

28 

co2 

0.0A2 

AA 

°2 

0.002 

32 

NO 

0.001 

30 

OH 

Trace 

17 

n2o 

Trace 

AA 

Calculated  composition  assuming  thermodynamic  equilibrium  at  the 
combustion  temperature.  Results  given  in  NASA  Handbook  and  quoted  by 
Pickett  et  al . '9^ 

TABLE  2 

Composition  of  the  Low-Earth  Orbit  Atmosphere ) 

(Molecules/cm^ ) 


Species 

200 

Altitude  (km) 
A00 

600 

H 

1. IxlO5 

0 

4.6xl09 

A. 8xl07 

8. 3x10s 

n2 

1. 7xl09 

6. 2x10® 

5xl02 

°2 

1.2x10® 

1.5xl05 

A.  A 

7 


Figure  1.  An  arbitrary  viewing  direction  repreacnted  by  a  unit  vector,  Tv> 
within  the  SOCRATES  cartesian  coordinate  system.  q  is  the  angle  between 
the  +X  direction  and  Tv,  and  <9  ia  the  angle  between  the  +Y  direction  end 
the  projection  oE  1  into  the  Y-Z  plane. 
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SCATTERING  RESULTS 


3  . 

Figures  2-3  contain  panels  of  the  H2O  column  density  for  the  nine 
calculations.  Figures  2  and  3  show  the  column  densities  obtained  by 
integrating  along  the  Y  direction,  in  gray  scale  and  contour  form, 
respectively.  Figures  4  and  5  show  the  column  densities  obtained  by 
integrating  along  various  directions  from  the  shuttle  bay,  as  a  function  of 
the  look  angles  a  and  0  defined  in  the  previous  section.  These  figures 
make  it  clear  that  the  choice  for  the  look  angles  provides  results  which 
are  independent  of  p  for  the  axisymmetric  cases.  In  the  spatial  plots 
(Figures  2  and  3),  the  atmospheric  wind  is  approaching  from  the  right.  The 
three  left  panels  show  the  altitude  variation  for  aligned  (0°)  plumes.  The 
effect  of  the  atmospheric  wind  increasingly  confining  the  plume  as  the 
altitude  is  decreased  can  be  seen  clearly  in  the  Figures,  even  with  the 
varying  length  scales.  For  the  90°  angle  of  attack  cases  (the  center 
columns)  the  plume  is  blown  back  substantially  within  the  first  kilometer 
at  200  km  altitude,  but  it  takes  about  ten  kilometers  for  substantial 
turning  to  be  evident  at  400  km.  At  600  km  the  plume  seems  largely 
unaffected  by  the  atmosphere.  Qualitatively  similar  results  are  evident 
for  the  retrofire  (180°)  cases  shown  in  the  right  column.  Figures  6-9  and 
10-13  show  the  corresponding  information  for  CO2  and  H2,  respectively. 

One  aspect  of  the  rarefied  nature  of  these  flows  is  that  light 
molecular  weight  species  tend  to  diffuse  much  more  readily  than  the  heavier 
species.  To  some  extent  this  is  evident  in  comparing  the  H2  to  the  H20  or 
C02  column  densities  in  the  previous  figures.  The  point  can  be  seen  more 
clearly  in  Figure  14,  which  compares  the  H2  and  H2O  number  densities  at  a 
position  2.1  km  downstream  from  the  exit  plane  for  the  200  km,  0°  case. 
Although  the  exhaust  contains  approximately  twice  as  much  H2O  as  H2  (Table 
1),  [H20j  -  2.8(82!  on  the  axis,  and  0.39[H2J  on  the  edge  of  the  plot.  The 
enhanced  diffusion  of  the  H2  causes  it  to  be  relatively  depleted  in  the 
dense  regions,  and  relatively  more  dominant  in  the  wing  regions.  As  will 
be  seen,  this  behavior  can  also  have  implications  for  radiation 


contamination. 
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Figure  2.  A  panel  of  gray  scale  plots  showing  the  H?Q  column  density 
resulting  from  VRCS  firings  at  200,  UOO,  and  600  Km,  each  at  three  angles 
of  attack:  0,  90,  and  I80u.  ~ 
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Figure  3.  A  panel  of  contour  plots  showing  the  H?0  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  Km,  each  at  three  angles 
ol  attack:  W,  90 ,  and  180°. 
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Figure  5.  A  panel  of  contour  plota  showing  the  H70  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  K»,  each  at  three  angles 
o(  attack:  0,  90,  and  180u.  The  column  density  ia  shown  as  a  function  of 
look  angles  for  different  viewing  directions  from  the  shuttle  bay. 
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Figure  6.  A  panel  of  gray  scale  plocs  showing  the  C07  column  density 
resulting  from  VRCS  firings  at  200,  600,  and  600  Km,  each  at  three  angles 
of  attack:  0,  90,  and  180^7  - 
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Figure  7.  A  panel  of  contour  plots  showing  the  CO?  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  Km,  each  at  three  anglea 
of  attack:  0,  90,  and  180u. 
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Figure  8.  A  panel  of  gray  scale  plots  showing  the  CO?  column  density 
resulting  from  VRCS  firings  at  200,  &00,  and  600  Ka,  each  at  three  angles 
of  attack:  0,  90,  and  180u.  The  column  density  ia  shown  as  a  function  of 
look  angles  for  different  viewing  directions  from  the  shuttle  bay. 
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Figure  9.  A  panel  of  contour  plot*  showing  the  CO?  column  density 
resulting  from  VRCS  firing*  at  200,  400,  and  600  Km,  each  at  three  angles 
of  attnck:  0,  90,  and  180^.  The  column  density  is  shown  ag  a  function  of 
look  angles  for  different  viewing  directions  from  the  shuttle  bay. 
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Figure  10.  A  panel  of  gray  scale  plots  showing  the  H?  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  Km,  each  at  three  angles 
of  attack:  0"j  9C ,  and  180^. 
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Figure  11.  A  panel  of  contour  plots  showing  the  H7  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  Km,  each  at  three  angles 
of  attack:  0,  90,  and  180u. 
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Figure  12.  A  panel  of  gray  scale  plots  showing  the  H?  column  density 
resulting  from  VRCS  firings  at  200,  400,  and  600  Km,  each  at  three  angles 
of  attack:  0,  90,  and  180u.  The  column  density  is  shown  as  a  function  of 
look  angles  for  different  viewing  directions  from  the  shuttle  bay. 
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Figure  13.  A  panel  of  contour  plota  showing  the  H7  column  density 
resulting  from  VRCS  firings  at  200,  4QQ(  and  600  Km,  each  at  three  angles 
of  attack:  0,  90,  and  180^.  The  column  density  is  shown  as  a  function  of 
look  angles  for  different  viewing  directions  from  the  shuttle  bay. 
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An  important  channel  for  absorbing  the  high  collision  energy  between 
rhe  ambient  0  atoms  and  the  plume  is  the  ccllisional  excitation  of 
vibrational  modes  of  plume  molecules  (particularly  the  heteronuclear 
molecules),  e.g. 

°f ast  +  »2°  -  H20(v3-1)  +  0fllow  (2) 

HjQCvj-D  -  H20  +  hv  (2.7  m)  (3) 

This  reaction  has  been  studied  both  theoretically  and  experimentally,  and 
found  to  be  very  fast^17-*^  In  this  work  we  have  used  the  value  k2  - 
[1. 2xlO_16T1-34e(_10730/kT) J  (Ref.  5).  Using  this  rate  coefficient  and  the 
known  17  msec  radiative  lifetime  of  H20(v3»l)  we  can  calculate  the  contours 
of  the  2.7m  radiation  caused  by  the  interaction  of  fast  0  atoms  with  H20 
from  the  plume.  This  is  shown  in  Figures  15-18,  again  for  all  the  cases 
considered.  In  a  similar  fashion  we  can  calculate  the  emission  in  the  4.3m 
band  excited  by  the  collision  between  fast  ambient  0  atoms  with  C02  in  the 
exhaust  of  the  VRCS  engines: 


°fast  +  C02 

-  C02(v3-1)  +  °8low 

(*> 

CO2  (  V3-I  )  *♦ 

C02  +  hv  (4.3m) 

(5) 

Here,  the  rate  coefficient,  k^,  is  less  well  known  and  we  have  used  k^  - 
(3. 5x10~16T2- 18e^  6.74/kT)j  (Ref.  5).  The  results  are  presented  in  Figures 
19-22. 

The  gray  scale  utilized  in  these  figures  encompasses  five  orders  of 
magnitude,  versus  two  for  the  column  density  figures  of  the  previous 
section.  The  reason  for  this  is  that  there  is  a  well  defined  vacuum  limit 
for  species  densities,  whereas  the  col 1 isional ly  induced  emission  depends 
on  an  interaction  with  the  atmosphere,  and  becomes  increasingly  tenuous  as 
the  altitude  is  increased.  The  spatial  extent  of  the  emission  at  0  and  90 
degrees  angle  of  attack  is  greater  than  the  density,  however,  even  without 
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Figure  16.  A  panel  of  contour  plots  shoving  the  H70  Vj  emission  at  2.7  um 
due  to  excitation  from  collisions  with  atomic  oxygen. 
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Figure  17.  A  panel  of  gray  scale  plots  showing  the  H70  emission  at  2.7 
pm  due  to  excitation  from  collioiona  with  atomic  oxygen.  The  emission  is 
shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  boy. 
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Figure  18.  A  panel  of  contour  plots  showing  the  H7Q  vj  emission  at  2,7  pm 
due  to  excitation  From  collisions  with  atomic  oxygen.  The  emission IT 


shown  as  a  function  of  look  angles  for  different  viewing  directions  from 


the  shuttle  bay. 
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Figure  19.  A  panel  of  gray  scale  plots  showing  the  C07  emission  at  4.3 
mdi  due  to  excitation  from  collisions  with  atomic  oxygen. 
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Figure  20.  A  panel  of  contour  plota  showing  the  CO?  Vj  emiasion  at  4.3  Mm 
due  to  excitation  from  collisiona  with  atomic  oxygen. 
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Figure  21.  A  panel  of  gray  scale  plots  showing  the  CO?  emission  at  4.3 
iia  due  to  excitation  from  collisions  with  atomic  oxygen.  The  emission  is 
shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 
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Figure  22.  A  panel  of  contour  plots  showing  the  CO?  emission  at  A.  3  mir 
due  to  excitation  from  collisions  with  atomic  oxygen.  The  emission  is 
shown  as  a  function  ot  look  angles  for  different  viewing  directions  from 
the  shutcle  bay. 


the  scale  change.  This  is  because  the  small  portion  of  the  flow  which  haa 
a  relatively  large  initial  upstream  velocity  component  suffers  higher 
energy  collisions,  and  is  more  likely  to  produce  emission.  The  importance 
of  this  portion  of  the  flow  is  enhanced  for  collisionally  induced  emission 
relative  to  density. 

The  total  solution  region  emission  for  these  processes  is  presented  in 
Table  3.  It  can  be  seen  that  the  effect  of  angle  of  attack  is  substantial 
for  both  mechanisms,  while  altitude  has  a  lesser  effect  between  200  and  400 
km.  For  the  600  km  cases,  there  is  substantial  radiation  outside  of  the 
solution  region,  though  it  is  spread  over  such  a  vast  region  that  the 
neglected  portion  is  not  important. 

Table  3 

Total  Solution  Region  Emission  (Watts/Sr) 


Case 

(Alt.  ,  A.  0.  A) 

H20(v3) 
(2.7  p) 

C02 ( v3 ) 
(4.3  m) 

200  Km, 

0° 

29. 

0.28 

200  Km, 

90° 

130. 

1.7 

200  Km, 

180° 

2  50. 

4.0 

400  Km, 

0° 

15. 

0.  14 

400  Km, 

90° 

160. 

2.3 

400  Km, 

180° 

310. 

6.2 

600  Km, 

0° 

0.76 

0.0082 

600  Km, 

90° 

14. 

0.22 

600  Km, 

180° 

29. 

0. 58 
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5.  REACTIVE  COLLISIONS 


We  consider  in  this  section  the  possibility  that  reactions  which  do 
not  proceed  at  room  temperature  due  to  a  large  activation  energy  might 
become  important  in  low  earth  orbit  because  of  the  orbital  velocity.  The 
reactions  which  we  consider  here  represent  species  which  are  important 
exhaust  products  of  VRCS  engines: 


°f  ast 

+  H20  -  OH*  ♦  OH 

(6) 

^*f  ast 

+  H2  -•  OH*  +  H 

(7) 

®f  ast 

+  co2  CO'  +  o2 

(8) 

There  are  some  high  temperature  measurements  for  all  three  reactions. 
However,  before  we  present  the  results  of  the  calculations,  a  discussion  of 
each  of  the  reactions  and  the  rate  coefficients  used  is  warranted. 

Reaction  (6):  The  reaction  to  produce  ground  state  OH  is  endothermic  by 
0.73  eV,  but  at  orbital  velocity  it  has  an  excess  energy  of  1.86  eV  (ECM 
for  0  +  HjO  at  7.7  km  s-*  is  2.59  eV  if  the  velocity  of  the  exhaust  is 
ignored.  If  the  velocities  are  added,  i.e.  at  e  -  180°,  then  the  reaction 
has  an  excess  energy  of  4.5  eV) .  Since  «e  for  OH  -  3738  cm"*  (2.7  m  or 
0.46  eV),  OH*  can  be  formed  at  v''  a  3  for  the  scenario  where  the  exhaust 
is  at  90°  with  respect  to  direction  of  travel.  At  180°  the  reaction  can 
excite  OH  into  the  A  (^E+)  state,  which  lies  at  32684  cm-*  (4.03  eV)  above 
the  ground  state.  (^0)  This  possibility  would  lead  to  emission  of  intense 
ultraviolet-visible  light,  since  the  A-»X  transition  is  allowed,  the 
radiative  lifetime  being  0.7  ms. 

The  rate  coefficient  for  this  reaction  has  not  been  measured  at 
energies  corresponding  to  orbital  energies.  However,  Albers  et  al.^*' 
have  reported  high  temperature  measurements  which  yield  a  rate  coefficient, 
k^,  of  -6.6  x  10“**e  (17300/kT)  cm3  molecule-*  a-*.  For  the  purpose  of 
this  work  we  have  taken  the  available  rate  constant  and  adjusted  the 
activation  energy  to  account  for  the  energy  of  the  first  vibrational  state 
and  (as  a  separate  reaction)  the  A  (^E+)  state. 
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Reaction  (7):  High  temperature  measurements  by  Miyauchi  et  al.^22^  give  k7 
-  3  y  1 n“ l^Te- (^480/T)  cm-3  molec  s-1.  Other  Measurements  have  been 
reported  by  Schott  et  al.,^22^  Pamidimukka la  and  Skinner , (2<4)  and  by 
Marshall  and  Fontijn^251  with  slightly  varying  results.  We  have  used  the 
data  of  Miyauchi  et  al.  The  reaction  is  endothermic  by  0.09  eV,  while  in 
low  earth  orbit  it  has  an  excess  energy  of  0.45  eV  at  90°  collision  and  0.9 
eV  at  180°  collisions.  Thus  only  vibrationally  excited  OH  can  be  formed. 
For  the  purpose  of  this  calculation  we  have  also  adjusted  the  activation 
energy  to  reflect  the  energy  of  the  first  vibrational  state. 

Reaction  (8):  This  reaction  is  endothermic  by  0.33  eV.  At  orbital  velocity 
and  at  90°  collision  angle  it  has  an  excess  energy  of  3.26  eV,  while  at 
180°  collision  angle  it  has  an  excess  energy  of  4  eV.  In  either  case,  the 
energy  is  not  enough  to  form  either  02  or  CO  in  excited  electronic  states. 
Vibrational  excitation  of  02  is  not  important  since  it  cannot  radiate  this 
energy.  This  means  that  vibrational  excitation  of  CO  is  the  only  channel 
left  for  important  internal  excitation.  Shock  tube  measurements  yield  a 
rate  coefficient  of  2  x  10-11  cm^  molecule-1  s-1  (Baber  and  Dean^26^).  We 
have  again  adjusted  this  rate  coefficient  to  take  into  account  the 
activation  energy. 

Table  4  presents  a  summary  of  the  rate  coefficients  used  in  these 
calculations  and  thermochemical  data  for  the  individual  reactions.  Figures 
23-26  show  the  results  for  reaction  (6a)  (with  vibrationally  excited  OH  as 
a  product),  which  can  be  compared  with  Figures  27-30  showing  the  results  of 
reaction  (7).  Note  that  in  these  figures  the  scale  was  changed  between  the 
two  reactions  to  reflect  the  lesser  emission  resulting  from  reaction  (7). 
This  can  be  seen  more  easily  in  Figure  31,  for  the  single  case  of  200  km 
altitude  and  0°  angle  of  attack.  The  lesser  amount  of  H2  in  the  plume  and 
the  higher  activation  energy  make  the  mechanism  of  reaction  (7)  generally 
less  effective.  The  exception  to  this  rule  is  when  wide  angle  areas  are 
considered  and  H2  may  be  much  more  effective  at  arriving  at  the  location 
prior  to  excitation. 

Figures  32-35  show  the  total  OH  emission  resulting  from  reactions  (6a) 
and  (7),  which  can  be  compared  with  the  result  of  the  generally  dominant 
H20  mechanism  in  Figures  23-26.  The  scales  are  the  same  in  these  two  sets 
of  figures,  and  it  can  be  seen  that  the  total  OH  emission  curves  show 
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Figure  23.  A  panel  of  gray  scale  plots  showing  the  OH  vibrational  emission 
at  2.7  Mm  as  a  result  of  the  reaction  of  0  +  H20  -*  OH*  +  OH.  ”™  *”  ™ 
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Figure  24.  A  panel  of  contour  plots  showing  the  OH  vibrational  emission  at 
2 . 7  pm  as  a  result  of  the  reaction  of  0  +  H20  ~  OH*  +  OH. 
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Figure  25.  A  panel  of  gray  scale  plots  showing  the  OH  vibrational  emission 
at  2.7  pm  as  a  result  of  the  reaction  of  0  +  H20  -»  OH*  +  OH.  The  emission 
is  shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 
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Figure  26.  A  panel  of  contour  plots  ahowing  the  OH  vibrational  emission  at 
2.7  pm  as  a  result  of  the  reaction  of  O  ♦  H20  -  OH*  +  OH.  The  emission  is 
shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 


-  38 


2  position  mm  z  position  mm  z  position  mm 


OH  EMISSION  AT  2.7  U*  FROM  0+H2 


600  KM,  0  DEGREES 


600  KM,  90  DEGREES 


600  Kn,  160  DEGREES 


-15-10  -5  0  5 

x  position  mm 

M00  Kfl,  0  DEGREES 


-20-10  0  10  20  30 

x  position  mm 

M00  KM,  90  DEGREES 


-50  -25  0  25  50 

x  position  mm 

>400  Ktl,  180  DEGREES 


-15-10  -5  0  5 

X  POSITION  <KM) 

200  «n,  0  DEGREES 


■  r 

C  l  ■  !  *!■  ^  •  • 


-S  -R  -3  -2  -1  0  1 

x  position  mm 


-20-10  0  10  20  30 

x  position  mm 

200  Ktl,  90  DEGREES 


-S0  -25  0  25  50 

x  position  mm 

200  Ktl,  100  DEGREES 


IV 

cn 

,  00* 

10' 

9 

=  i. 

.39* 

10' 

■9 

=  3. 

.  87x 

10' 

■10 

=  1 . 

,08x 

10' 

■10 

=  3. 

.  00X 

10 

•11 

=  8. 

,3Mx 

10' 

■12 

-2. 

.32x 

10 

•12 

=  6 

.  46x 

10 

-13 

-1 

.80x 

10 

■13 

<5 

,00x 

10 

-It 

WRTTS/SR/CM‘ 


-5-0-3-2-10  I  2 

x  position  mm 


x  position  mm 


Figure  27.  A  panel  of  gray  scale  plots  showing  the  OH  vibrational  emission 
at  2.7  pm  as  a  result  of  the  reaction  of  0  +  H2  -»  OH*  +  H. 
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Figure  28.  A  panel  of  contour  plots  showing  the  OH  vibrational  emission  at 
2.7  mw  a»  a  result  of  the  reaction  of  0  +  H2  —  OH”  +  H. 
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Figure  29.  A  panel  of  gray  scale  plots  showing  the  OH  vibrational  emission 
at  2.7  mm  as  a  result  of  the  reaction  of  0  +  H2  -»  OH*  +■  H.  The  emission  is 
shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 
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Figure  30.  A  panel  of  contour  plots  showing  the  OH  vibrational  emission  at 
2.7  um  as  a  result  of  the  reaction  of  0  +  H2  -»  OH* " ♦  H.  The  emission  is 


shown  as  a  function  of  look  angles  for  different  viewing  directions  from 


the  shuttle  bay. 
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Figure  32.  A  panel  of  gray  scale  plots  showing  the  total  OH  vibrational 
emission  at  2. 7  mb  as  a  result  of  the  reactions  O  +  H20  -»  OH*  »  OH  and  0  + 
H2  -  OH*  +  H. 
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Figure  33.  A  panel  of  contour  plots  showing  the  total  OH  vibrational 
emission  at  2.7  pm  as  a  result  of  the  reactions  0  H20  -•  OH*  +  OH  and  0 
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Figure  34.  A  panel  of  gray  scale  plots  showing  the  total  OH  vibrational 
emission  ac  2.7  Mm  as  a  result  of  the  reactions  0  +  H20  -»  OH*  +■  OH  and  O  + 
H2  -»  OH*  +  H.  The  emission  is  shown  as  a  function  of  look  angles  for 
different  viewing  directions  from  the  shuttle  bay. 
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Figure  35.  A  panel  of  contour  plots  showing  the  total  OH  vibrational 
emission  at  2,7  nm  as  a  result  of  the  reactions  0  +  H20  -»  OH*  *  OH  and  0  + 
H2  -•  OH’11  +  H.  The  emission  is  shown  as  a  function  of  look  angles  for 
different  viewing  directions  from  the  shuttle  bay. 
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filled  in  regions  at  wide  angles  which  are  not  present  for  the  HjO 
mechanism  alone.  The  difference  between  the  two  curves  is  the  H2  mechanism 
which,  due  to  the  greatly  enhanced  diffusivity  of  H2,  shows  up  at  the  wide 
angles. 

In  comparing  Figures  15-18  with  Figures  32-35,  it  can  be  seen  that 
comparable  emission  is  predicted  for  the  reactive  mechanisms  (Reactions  6a 
and  7)  and  T-»V  mechanism  (Reaction  2).  This  suggests  that  the  reactive 
mechanism  should  be  considered  in  future  2.7p  emission  calculations. 

The  predicted  ultraviolet  emission  from  reaction  (6b)  (with 
electronically  excited  OH  as  a  product)  is  shown  in  Figures  36-39  for  all 
cases.  The  very  high  activation  energy  for  this  mechanism  results  in  an 
emission  profile  which  is  quite  different  from  previous  density  and 
emission  profiles.  In  particular,  the  200  km  0°  case  only  has 
contributions  from  the  very  small  portion  of  the  flow  which  is  predicted  to 
turn  and  head  back  upstream.  The  model  was  not  designed  to  be  accurate  for 
that  portion  of  the  flow,  so  such  results  should  be  considered  preliminary. 

The  predicted  CO  vibrational  emission  from  reaction  (8)  is  shown  in 
Figures  AO-43  for  all  cases.  In  comparing  these  results  to  the  direct  C02 
excitation  shown  in  Figures  19-22,  it  is  once  again  seen  that  the 
signatures  are  generally  predicted  to  be  comparable,  although  in  this  case 
the  wavelength  is  shifted  slightly  for  the  CO  vibrational  emission. 


TABLE  4 

Rate  Coefficients  Used  in  SOCRATES 
k  =  ATnexp( -Ea/kT) 


REACTION 
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Ea 

(kcal/mole) 

AH 
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H2O  -*  0  +  H2CK  Vj ) 

1 . 2 (-16) 
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CO2  “*  0  ♦  CO2  (  Vj  ) 

3 . 5 (-2 1 ) 

2.18 
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110. 20 

109.  7 

7. 

0 

•f 
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3. 0(-14) 
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14.30 
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Figure  36.  A  panel  of  gray  scale  plota  showing  the  OH  ultraviolet  emission 
at  306.4  nm  as  a  result  ot  the  reaction  0  +  H?0  ->  OH(A)  +  OH. 
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Figure  37.  A  panel  of  contour  piota  showing  the  OH  ultraviolet  emission  at 
306.4  run  as  a  result  of  the  reaction  0  +  H70  -»  OH(A)  ♦  OH. 
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Figure  38.  A  panel  of  gray  scale  plots  showing  the  OH  ultraviolet  emission 
at  306.4  nm  as  a  result  of  the  reaction  0  +  H70  -»  OH(A)  »  OH.  The  emission 
is  shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 
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Figure  39.  A  p.inel  of  contour  plots  showing  the  OH  ultraviolet  emission  at 

306.4  nm  as  a  result  of  the  reaction  0  ♦  HjO  -»  OH(A)  ♦  OH.  The  emission  is 

shown  as  a  function  of  look  angles  for  different  viewing  directions  from  « 

the  shuttle  bay. 
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Figure  40.  A  panel  of  gray  scale  plots  showing  the  CO  vibrational  emission 
at  4.6  m<o  as  a  result  of  the  reaction  of  O  ♦  C02  -»  CO*  +  02. 
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Figure  41.  A  panel  of  contour  p!ot«  showing  the  CO  vibrational  emission  at 
4.6  mw  aa  a  result  of  the  reaction  of  0  +  CQ2  -*  CO*  +  02. 
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Figure  42.  A  panel  of  gray  scale  plota  showing  the  CO  vibrational  emission 
at  4.6  mw  aa  a  result  of  the  reaction  of  0  +  CQ2  —  CO*  +  02.  The  emission 
is  shown  as  a  function  of  look  angles  for  different  viewing  directions  from 
the  shuttle  bay. 
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Figure  43.  A  pnnel  of  contour  plots  ahowing  the  CO  vibrational  emission  at 
4,6  urn  as  a  result  of  the  reaction  ot  0  +  CQ2  -  CO*  +  02.  The  emission  ia* 
shown  as  a  function  of  look  angles  for  different  viewing  directions  trow 
the  shuttle  bay. 
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6. 


SUMMARY 


SOCRATES  is  already  an  effective  tool,  which  will  become  even  more  so 
as  the  planned  additional  capabilities  are  added.  It  has  been  used  to 
systematically  investigate  the  effects  of  altitude  and  angle  of  attack  on  a 
variety  of  scattering,  excitation  and  reactive  mechanisms.  A  major  new 
result  from  this  work  is  that  reactive  mechanisms  are  predicted  to  be  aa 
important  as  collisional  excitation  in  producing  IR  radiation,  and 
potentially  important  in  the  UV,  as  well.  The  approach  outlined  here  makes 
it  possible  to  contemplate  making  optical  meaaurements  from  the  Space 
Shuttle  or  other  large  space  platforms  and  planning  to  identify  and 
subtract  contributions  from  the  local  environment  to  the  observed  signals. 
A  major  requirement  is  that  some  information  be  available  about  the 
composition  of  the  gaseous  local  environment. 
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